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Abstract. In this paper we describe invariant geometrical structures in the phase space of the Swift-Hohenberg 
equation in a neighborhood of its periodic stationary states. We show that in spite of the fact that these states are 
only marginally stable {i.e., the linearized problem about these states has continuous spectrum extending all the way 
up to zero), there exist finite dimensional invariant manifolds in the phase space of this equation which determine the 
long-time behavior of solutions near these stationary solutions. In particular, using this point of view, we obtain a new 
demonstration of Schneider's recent proof that these states are nonlinearly stable. 
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1. Introduction 

In this paper, we study the non-linear stability of space-periodic, time-independent solutions of 
the Swift- Hohenberg equation 

d t u = {e 2 -(l + d 2 x ) 2 )u-u\ (1.1) 

Here, u(x,t) is defined on R x R + and takes real values and e > is a small parameter. The 
Eq.(l.l) has stationary solutions u(x, t) = it (x) which are of the form 

n6Z 

The non-linear stability problem addresses the question of the time evolution of initial data 
which are close to u e w , and stability in this context means that the solution converges to u £ u 
as t — > oo. The range of possible values of u; is given by e 2 > (1 — u 2 ) 2 when cj is close to 
1. To simplify the exposition we shall concentrate on the case to = 1, and omit henceforth the 
index to. 

In a very interesting paper, G. Schneider [Sch] has solved this problem, and the present 
work relies heavily on his ideas. Our aim is to simplify somewhat the exposition of [Sch] and 
to extend the result by giving a more precise asymptotic analysis, using the description of the 
asymptotic behavior in terms of a continuous renormalization group and invariant manifolds as 
introduced in [W], see below. 

The existence of solutions of the form Eq.(1.2) is a well-established fact, (see e.g. [CE]) 
and we repeat here only those points of the discussion which are needed in the sequel. The 
equation for the stationary solution is F(it, e) = 0, where 

F(u, e) = (e 2 - (1 + d 2 x ) 2 )u - u 3 . (1.3) 

The equation F = has the trivial solution it = 0, e = 0. Linearizing around this solution, we 
see that DF equals 

DF = -(1 +<9 2 ) 2 ©0, 
acting on some weighted subspace of F 2 (R) © R. The null space of DF is spanned by 

{cosx, sinx} © and 0©1, (1.4) 

and thus, bifurcation theory suggests the existence of solutions of the form of Eq.(1.2), when 
e 7^ 0. This is indeed what happens (cf. [CR], [CE]), and the higher frequency terms in Eq.(1.2) 
are generated from the basis Eq.(1.4) by the non-linearity u . The method clearly extends to 
similar polynomial non-linearities. An explicit calculation shows that F(it e , e) = for 



2 

u s (x) = e— =cos(x)+£ h £ (x) , 



(1.5) 
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and h £ (x) = h s (x + 2tv). Thus, the function u £ equals u £ l of Eq.(1.2). We have broken the 
translation invariance of the problem by the choice of cos in Eq.(1.5), instead of, say, sin. 

We next pass to the linear stability analysis of the solution u £ . This is again a classical sub- 
ject, initiated by Eckhaus [E], which we summarize for convenience, see also [CE]. Linearizing 
Eq.(l.l) around the solution u £ we are led to study the operator L £ = (e 2 — (1 + <9 2 ) 2 ) — 3u 2 £ , 
that is, 

(L £ v)(x) = (e 2 - 3u 2 £ (x))v(x) - (1 + d 2 x ) 2 v(x) . 

Because u £ is a In periodic function, it is most convenient to work in Floquet coordinates 
{i.e., with Bloch waves). To fix the notation, we give some details: Begin by introducing the 
following representation for / G L 2 (R): 

r ^ r 1 / 2 

fix) = J dke~ lkx f(k) = J dee- imx e- Ux f(m + e) 



y l 2 



die-^f e (x) 

1/2 



where 

f e (x) = J2 e ~ imX f(m + £). (1-6) 

Properties of /. Note first that f e is 2n periodic. Furthermore, the definition of f e (x) can 
be extended to all £ G R by the definition 

f e+ i(x) = e~ lx f e (x) . 

We next observe that if / has a smooth, rapidly decaying Fourier transform, then f e (x) will also 
be a smooth function of £ and x. If /, g are in L 2 (R), then it follows from the definition of f t 
that ^ 

(fg)7(x) = [ dkf e _ k (x)~g k (x). (1.7) 

7-1/2 

We finally note that if s is a 2tt periodic function, then 

g t (x) = 8(£)s(x) . (1.8) 

It is now easy to see that 

{L e v)~{x) = (e 2 -(l + (t£ + d x ) 2 ) 2 )v e (x)-3(u 2 £ v)7(x). 

In the language of condensed matter physics, £ is the quasi-momentum in the "Brillouin zone" 
[— i, and L £ leaves the subspace J g of functions with quasi-momentum £ invariant. Using 
the properties just described, we get 

= {e 2 -(l + (z£ + d x ) 2 ) 2 )v e (x)-3u 2 £ (x)-v e (^ = {L s , e v £ )(x) . (1.9) 
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To fix the notation, we repeat the calculation done by Eckhaus, cf. also [CE], [M]. We denote 
c(x) = cos(x), s(x) = sin(x). The method of Eckhaus consists in projecting the eigenvalue 
problem for L e g onto the subspace spanned by the "bifurcating directions" c and s. Observe 

that, modulo higher frequency terms, we have c 3 = |c, c 2 s = |s, and therefore the projection 
of L el onto this subspace is described by the matrix 

f-4£ 2 -£ 4 -2e 2 + 0(e 4 ) -Ai£ 3 \ , 4) f 0{£ 2 ) G(£) \ 
^ Ai£ 3 -4£ 2 -t) +U ^ £ > { 0(£) 0(£ 2 ) J ■ 

The eigenvalues of this matrix are 

A°,o = -(4 + 0(e 2 ))£ 2 + 0(£ 3 ), 

A° t = -2(e 2 + 0(e 4 )) - (4 + 0(e 2 ))£ 2 + 0(£ 3 ) + 0(£ 4 + e 4 ) . 

Thus, the restriction of L e e on the subspace spanned by c and s has its spectrum in the left 
half -plane. Note that the corresponding eigenvectors are s + 0(£ + e) and c + 0(£ + e). 
Extending this calculation to the full space, one shows in the same way [E, CE, M] that 

Theorem 1.1. For sufficiently small e > the operators L e e , with£ E [— i, ^] are selfadjoint 
on the Sobolev space H 4 , have compact resolvent and a spectrum satisfying 

\o(e 2 ) = -(4 + 0(e 2 ))£ 2 + 0(£ 3 ) = -c Q (e 2 )£ 2 + 0(£ 3 ) , 
X tA (e 2 ) = -2{e 2 + 0(e 4 ))-{4 + 0(e 2 ))£ 2 + 0(£ 3 ), (1.10) 
X ej < -(l-/) 2 + 0(£ 2 ), j =2,3,... . 

Notation. Since we mostly concentrate on the branch 0, we shall abbreviate X e = X e (e 2 ). 
The eigenfunction corresponding to X e is 

ip e e (x) = const. ( u' e (x) + i£g s (x) + h s ^(x)£ 2 ) , 

where u s is the stationary solution, and both g s and h e £ are 27r periodic. If we choose the 

constant to normalize the L 2 norm of <p s e to 1, then <p £ £ = ix~ 1 / 2 sin(x) + 0{e + \£\). 

We can now formulate the main question of this paper: Having seen that the solution u e is 
linearly (marginally) stable, is it true that this solution is stable under the non-linear evolution? 
The answer will be affirmative. As pointed out by Schneider [Sch], the result is not obvious, 
since the leading non-linear term does not have a sign. Indeed, the non-linear evolution equation 
for a (small) perturbation of u e is 

d t v = -(1 + d 2 x ) 2 v + e 2 v - 3u 2 e v - 3u s v 2 - v 3 , 

where we recall that u s is of order e, and approximately equal to O(e) cos(x). Reducing again 
to quasi-momentum £, and using Eq.(1.8), we get the equation 



(1.12) 



Nonlinear Stability 



5 



and it is the term 3u s (v 2 )^ which does not have a sign. The saving grace will be the diffusive 
behavior suggested by the spectrum (in particular the branch X e ). At first sight, the non-linearities 
seem to be too singular for diffusion to dominate a potential divergence. Indeed, it is well known 
that, e.g., the equation 

d t u = d 2 x u + u 3 , 

has solutions which blow up in finite time [L], and the quadratic term makes things even worse. 
The beautiful observation of Schneider[Sch] is, however, that the problem Eq.(1.12) is rather of 
a form reminiscent of 

d t v = d 2 x v - d 2 x (v 2 + v 3 ) , (1.13) 

which is good enough for convergence [CEE, BK, BKL]. 

In later sections we examine in detail the form of the non-linear terms in Eq.(1.12), but 
here we explain briefly why these terms are similar to the non-linear terms in Eq.(1.13). The 
derivatives in the non-linearity have their origin in the symmetries of the problem, and they 
are easier to understand in momentum space. In fact, Eq.(1.13) is a good approximation to 
Eq.(1.12) only in the low-momentum (small £) regime, but this is sufficient since for £ outside 
a neighborhood of £ = 0, the stationary solutions are linearly stable, (and not only marginally 
stable) and the form of the non-linearity is unimportant. 

To understand the low-momentum behavior of Eq.(l . 1), note first that the Swift-Hohenberg 
equation Eq.(l.l) — and, incidentally, other equations with coordinate independent right hand 
side — has a circle of fixed points generated by translations. If we now study the Eq.(1.12) at 
£ = 0, this corresponds to studying the Swift-Hohenberg equation in the space of functions 
of period 2n. In this space, say L 2 ([0, 2-7r]), the linear operator in Eq.(1.12) has pure point 
spectrum with a simple eigenvalue at and all other eigenvalues real and strictly negative. In 
this case, as Schneider notes, the center manifold theorem can be applied, and there exists a 
1 -dimensional center manifold. We also see immediately that the eigenvector corresponding to 
the eigenvalue is d x u e , i.e., it is tangent to the circle of fixed points generated by translations. 
In fact, since any fixed point sufficiently close to the origin must lie in the center manifold, we 
see that the center manifold coincides with the 1 -dimensional circle of fixed points. Thus the 
non-linearity in the equation, when restricted to the center manifold, must vanish. This shows 
that the effective non-linearity in Eq.(1.12), when evaluated at £ = 0, must vanish and this 
accounts for one derivative in Eq.(1.13). More precisely, we see that the effective non-linearity 
inEq.(1.12) is bounded by 0(£), as is the non-linearity inEq.(1.13). The second derivative of the 
non-linearity in Eq.(1.13) arises because of "momentum conservation." Since ip e e is a smooth 
function of £, the linear term in Eq.(l.l 1) must of the form i£g e , with g e independent of £. Since 
the interaction is local in x, one sees upon working out the integrals that all terms proportional 
to £ in the non-linearity cancel exactly, see Eq.(A.3). Thus, the low momentum behavior of 
Eq.(1.12) is as if the non-linearity was differentiated twice — i.e., exactly as in Eq.(1.13). 

Our main result is that this intuitive argument correctly predicts that the leading order 
asymptotics are diffusive, and that furthermore, the higher order asymptotics are controlled by 
a sequence of finite dimensional invariant manifolds. Thus, our approach provides some insight 
into how finite dimensional geometrical structures can arise from a problem with continuous 
spectrum. 
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Stability Theorem 1.2. Fix n > 1 and 5 > 0. There exists a Hilbert space, 7i(n), such 
that if the initial conditions of (1.12) lie in a sufficiently small neighborhood of the origin in 
H(n), then there exists an n + 1 dimensional, invariant manifold in the extended phase space 
P{n) = R+ x H(n) of (1.12), and any sufficiently small solution of this equation which is not 
on this manifold approaches it at a rate 0(t~^ n+1 ~ s ^ 2 ). In particular, ifn = 1, small solutions 
of (1.12) have the asymptotic form: 

Remark. In Sections 2 and 3, we will make clear precisely what the Hilbert spaces H(n) are 
and what we mean by "sufficiently small." 

The remainder of the paper is devoted to a proof of the Stability Theorem 1.2. 



2. Formulating the Stability Theorem 1.2 in terms of scaling vari- 
ables 

In this section, we transform the problem to a rescaled dynamical system. In the next section, 
we will cast the dynamical system thus obtained into an invariant manifold problem. 

The idea of the proof is to focus on the "central branch" of the spectrum, X £ = X £ (e 2 ), 
which is only marginally stable. The relevant part of the spectrum for the long-time asymptotics 
is only the part in a small neighborhood of I = 0, a fact we exhibit by an appropriate rescaling of 
the dependent and independent variables. This rescaling has the disadvantage that it introduces 
a singular perturbation in the variables corresponding to the "stable branches" of the spectrum, 
X en (e 2 ), n > 1, because the corresponding modes decay extremely fast, when rescaled (at 
least on a linear level). However, invariant manifold theory has long been used to treat singular 
perturbation problems, and we are able to use it for that purpose here as well. In addition, these 
invariant manifolds will provide us with a geometric description of the long-time asymptotics 
of solutions near the stationary states. 

Our method generalizes to other problems of similar spectral nature, see the example of a 
cylindrical domain given in [W2] . 

Henceforth, we fix e > 0, and omit it from most subscripts. Since L e = L £ e is self-adjoint, 

we can define the (orthogonal) spectral projections P £ and P^~, which project onto the central 
branch and its complement. 

Remark. We know that for \£\ sufficiently small, say \£\ < £ /2, one has 

spec(iy^) = -c (e 2 )f + 0(£ 3 ) , 

and that this is the eigenvalue closest to in spec(L^). We continue this projection smoothly to 
larger I even if it cannot be guaranteed to be a projection onto the highest eigenvalue. But note 
that for those values of £ the spectrum of L t can be shown to be strictly bounded away from 0, 
see, e.g., [CE, page 102]. 
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To study the non-linearity, and to show the mechanism leading to the result which is 
analogous to Eq.(1.13), we write the Eq.(1.12) in more detail: 



1 



,1/2 

d t M x ) = (L t v t )(x)-3u e (x) / dkv £ _ k (x)v k (x) 

J -1/2 

■V2 (2.1) 

dk x dk 2 v £ _ ki _ k2 (x)v k (x)v k2 (x) 

'-1/2 

= {L e v e )(x)-{F 2 (v)) e (x)-{F 3 (v)) £ (x). 

We now decompose the Eq.(2.1) by projecting onto P £ and P £ . If / G L 2 , we let f £ = P £ f £ , 
and f £ = P £ f £ . Similarly, L\ = P e L e P e and L £ = P^L £ P £ ^. Then we get 

d t v c £ {x) = m(x) - {P e F 2 (v) e ) (x) - {P e F 3 (v) e ) (x) , (2.2) 

and a similar equation for v^~: 

d t vl(x) = L\vi{x) - {P^F 2 (v) e )(x) - (PtF,{v)e){x) . (2.3) 

We next split the first equation into a piece corresponding to small i.e., \i\ < £ and 
another corresponding to large £. Since we want to construct invariant manifolds, we need some 
smoothness in this construction and we choose a smooth cutoff x satisfying 



1, if\£\<£ , 
0, if \£\ > 2£ , 



and of course £ < j. In fact, we shall choose £ > so small that P e is the projection onto the 
central eigenspace for all £ G [— £ , £ ]. Let (f £ denote the normalized eigenvector which spans 
the range of P e (for \£\ < £ , and smoothly continued for £ beyond that value). Then v\ can 
be written as v\ = V(£)cp £ , where it is understood that V is really a function of v. We also let 
U £ denote the operation Hgfi = (ip e \f e ), where (•) is the scalar product in 3 r £ . This operation 
extracts the coefficient V and therefore Eq.(2.2) can be written as 

d t V{£) = \ t V{£) - U £ P £ F 2 (v) e - U £ P £ F 3 (v) £ . (2.4) 

Defining V < (£) = X (£)V (£) , and V> (£) = (1 - X (t))V(£), the Eq.(2.4) can be rewritten as 

d t V<{£) = X e V<(£)-(r(V<,V>,v^))(£), 



where 



d t V>(£) = \ t V>{£)-(nv<,V>^)){£), (2 ' 5) 
{f(V<,V > ,v ± ))(£) = (l- X (£))(u e P e F 2 (v) e + U e P e F 3 (v) i 
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and 

v e (x) = (V<(£) + V>{£)) ■ <p t {x) + vi{x) . 

Note that since V > is supported outside [—£ , £ Q ], both it and v 1 ' decay exponentially (at least 
at the linear level) and hence will be irrelevant for the asymptotics of V < , as we shall show. 
With this in mind, we introduce a new coordinate, V s , which combines the "irrelevant" pieces, 
V s = (V > , -D -1 ). Then the Eq.(2.5) combined with Eq.(2.3) takes the more suggestive form 

d t V<(£) = X e V<(£)-(f(V<,V s ))(£), 
d t V s = L^V + giV^V), 

and we know that the spectrum of the linear operator is contained in ( — 00, — cr s ), for some 
a 8 > 0. 

In order to proceed further, we analyze the non-linear terms in Eq.(2.6) in more detail. In 
particular, we concentrate on the most critical terms, namely those in / of Eq.(2.6) which depend 

only on V<. We decompose f(V<, V s ) = / 2 ( ° ) (F < ) + / 3 ( ° ) (^ < ) + / 4 ( ° ) (^ < , ^ s ), where / 2 (0) 
collects the terms which are homogeneous of degree 2 in V < and those of degree 3. One 
gets 



{fl 0) (V<))(£) = 3 X (£) J dxJp e (x)u s (x) Jjk^ k {x)^_ k {x)V<{k)V<{£-k) 

r l/2 

= 3 X {£) / dkK 2 (£, k)V < (k)V < (£ - k) , 

J-l/2 



V2 (2.7) 



{ff\v<)){£) = X (£) [ dxip e (x) [ dk l dk 2 <p hl {x)<p hl {x)<p e _ hl _ hl {x) 

J 7-1/2 



1/2 

xV<(k l )V<(k 2 )V<(£-k 1 -k 2 ) 

rl/2 

= x {£) / dk l dk 2 K 3 (£,k l ,k 2 )V < (k l )V < (k 2 )V < (£- k x - k 2 ) . 
J -1/2 

At this point, we make use of the diffusive nature of the problem for V K , by introducing scaling 
variables as in [W]. This will give us a more precise description of the convergence process than 
the one obtained in [Sch]. We rescale the variables in Eq.(2.6) as follows: We first fix, once and 
for all, a (large) constant t > 0. Then we define 

V<(£,t) = w c (sign(e)J\A e \(t + t ),\og(t + t )) , 

1 , (2 - 8) 

V\l,t) = w s (sign(£)y/\A e \(t + t ),\og(t + t ))/(t + t ) 1 / 2 , 

where A e = X e for \£\ < £ /2 and is monotonically extended beyond that region in such a way 
that it is parabolic for large \£\. (This artifact is needed because we have no guarantee that X e 
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itself is monotone.) Note that if X £ were equal to —const, i , this scaling would amount to the 
usual "diffusive" rescaling. Our choice takes into account higher order corrections produced by 
higher order terms in X e . If we let now p = sign(£) ^\A £ \(t + t ), and r = log(t + t ), then 
Eq.(2.6) implies that w c and w s obey the following equations: 

d T w c = {-p 2 - \pd p )w c 

+ e T (/ 2 « e-/ 2 ) + f 3 (w c , e-/ 2 ) + f,{w\ w"e~^\ e"^ 2 )) , (2.9) 
e~ T d T w s = M exp( _ T/2) w s + \e~ T w* - \e~ T pd p w* + e r/2 g(w c , w s e~ T/ \ e" r/2 ) , 

where f 2 , / 3 , / 4 and M in Eq.(2.9) are defined below. If 



pe-/ 2 = p(t + t Q )- 1 / 2 = sign^y^, 
and if we denote the inverse transformation by 

£ = $(pe" r/2 ) , 



where $ is the inverse function of x t— > sign(a;)-\/jAJ, then, given a function iu = t), we 
define the nonlinearity 

[f 2 (w,e-r/ 2 )](p) = [fl°\w{-,e r ))](*(pe- T ' 2 )) 

= [/ 2 (0) («;(-,t + t ))]($(p(t + t )- 1 / 2 )). 

(Note that <&(x) = x(l + O(x)).) Analogous definitions apply to / 3 and / 4 . The operator M 
will be described in detail in Eq.(2.13). 

Remark. The non-linearities f 2 ,. . . depend on the choice of t . If we consider the initial value 
problem for the Swift-Hohenberg equation, the "smallness" assumption on the perturbation of 
the periodic state is to be understood with respect to a choice of a (sufficiently large) £ . As we 
will see, however, the nonlinear terms can be bounded, independent of t , for all t Q > T > 0. 

To this change of variables will correspond the following (non-exhaustive) list of substi- 
tutions in the integrals in Eq.(2.7): Let a, b E Then 

b r e T / 2 <^- l (b) 



po pe < v \o) 

X (£) dk -> X {®{pe- T/2 ))e- T ' 2 dq& \qe~^ 2 ) , 



K(M) -> «;(p,r), 

- w(A(p,g,r)). 



e T/2$-i(a) 

(2.10) 
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Here, we define 

T(p, q ,r) = $(pe" T / 2 ) - <%e" T / 2 ) , 

A(p, ? , r) = e^" 1 ($(pe" r/2 ) - $(ge" r/2 )) . 

It follows at once from the definition of $ that 

r(p, g, r) = e- r / 2 (;p - g) • (l + 7 (p, g, r)) , ^ 

A(p, g, r) = (p - g) • (1 + x(p, g, r)) , 

where x and 7 are bounded and smooth. 

We next discuss in detail the spectrum of M exp ^_ T ^ 2 y which is just the rescaled linear 

operator for the "stable" part of w, cf. Eq.(2.6). Recall first that V s = (V > ,v~ L ). This 
introduces a natural decomposition of tu s = (tu^tu^aswellasofM^^^ = M sxp ^_ T ^ 1© 
-^ex P (-r/2) 2- From the definition of the first component, we get 

(K q ,(-r/2),i/f)(P,r) = ^-(l + Cx+^Cpe-^ 2 )) 2 ) 2 -^^--/ 2 )))/^^), (2.13) 
where DC(£) is a kernel given by 

X(£) =3 dxlp e (x)u 2 e (x)ip e (x) . 



(Recall that (p £ really depends on e as well and should be written (p £ e .) Since V s has support 
bounded away from £ = 0, say \£\ > £ /2, we see that w\(p, r) will have support in \p\e~ T ^ 2 > 
^\A io / 3 \, and the spectrum of M exp ^_ T ^ 2 ) y i is seen to be contained in {cr|Recr < a < 0}, for 
some a Q and for all r > 0. 

A very similar argument detailed in Appendix B shows that the spectrum of -^ e xp(-r/2) 2 
is also contained in such a set. Thus, the linear evolution generated by A^ e xp(-r/2) contracts 
exponentially. See Lemma B.6 below for details. 

We next consider the operator L = (—p 2 — \pd p ), which appears in the first component 
of Eq.(2.9). The detailed study of the semi-group generated by L will be given in Appendix B. 
Here, we discuss its properties on an informal level. The Fourier transform of L is d 2 + \xd x + \ , 
which is conjugate to the harmonic oscillator H = d x — x /16 + 1/4 by the (unbounded!) 
transformation T, of multiplication by exp(x 2 /8). In formulas: L = T~ l H T. Therefore, H 
has (say, on L 2 ), discrete spectrum \i ■ = —j/2, j = 0, 1, . . . . It is this spectrum which leads 
to a nice interpretation of the convergence properties of the Swift-Hohenberg equation. The 
eigenvalues of L are unchanged by the transformation T, (and the eigenfunctions are multiplied 
by a Gaussian), so to each eigenvalue \i of L there corresponds a decay rate e TM in the linear 
problem. Because of the transformation of variables from t to r, this decay rate becomes 
(t + t ) M in the original problem Eq.(2.6). In other words: Neglecting the non-linearities in 
Eq.(2.9) and setting w s = 0, (and ignoring potential problems related to the unbounded operator 
T) we have a solution 

00 

w<( P ,t) = J2v- Tm/2 M m (2p), (2.14) 

m=0 
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where < K m is the m eigenf unction of L. In the original variables, this means that 

oo 

V<(£,t) = ^^ m (t + t )- m/2 ^ m (2^ + g 1 / 2 (l+O(|£| 1 / 2 ))) . (2.15) 

m=0 

Thus, to each m there corresponds a specific rate (fi m = —m/2) of decay for a part of the 
function V < . Note that a change of t just corresponds to a rearrangement of the series. (This is 
not contradictory, since a change of t also changes the initial condition, and hence the solution 
whose asymptotics we are computing.) In particular, the slowest rate of decay is associated 
with !K , which is Gaussian, and thus, at least at the linear level, a "generic" perturbation of the 
stationary state will decay like exp(— c£ 2 t), for some c > 0. In terms of the original independent 
variables (x, t), it decays like t~ x l 2 exp(— x 2 / {Ate)), as t — > oo. This means that at this level, 
the periodic stationary states are stable, and that perturbations of them decay like solutions of 
the linear heat equation. The invariant manifold theory guarantees that this behavior persists 
in the non-linear problem, and in fact it tells us more. We will see that in suitable spaces we 
can construct a sequence of manifolds M ■ of dimension j = 1,2,..., such that any solution of 
Eq.(2.9) approaches a solution on M • at a rate e TMj ~' , or again reverting to the original (x, t) 

variables, at a rate O((t+t )^~ l ). In the case at hand, this is O((t+t )~^ 2 ). Thus, in principle, 
we can analyze finer and finer details of the asymptotics of perturbations of the stationary state 
by considering the behavior of the solution on these finite dimensional manifolds. 

3. Casting the Stability Theorem 1.2 into an invariant manifold 
theorem 

At the end of the preceding section, we have seen that the spectrum of the linear part of Eq.(2.9) 
has the following nature: The component w c satisfies a differential equation whose linear part 
has eigenvalues ju • = —j/2,j = 0, 1, ... , N, provided we work on a space of sufficiently 
smooth and rapidly decaying functions. The evolution of w s is governed by an equation with an 
even more stable spectrum. 

The invariant manifold theorem will show in which sense the built-in scalings of Eq.(2. 14) 
survive the addition of non-linearities. While this presents no conceptual problems at all — and 
this is the beauty of the present approach — some care is of course needed in the application of 
the invariant manifold theorem. Another point which might be overlooked is the following: The 
invariant manifold theorem does not say that the representation of the full non-linear problem 
is the same as in Eq.(2.15), but with slowly varying w-. Rather, we will show that on the 
complement of a dimension j — 1 surface in the function space, the solutions decay at least like 
t~^ 2 , (for every j > 1), provided the initial data are sufficiently small and smooth. 

In order to apply the invariant manifold method to the problem, we need bounds on the 
non-linearities and bounds on the semi-group generated by L. While the factor of t = exp(r) 
in front of f 2 in Eq.(2.9) might look like a disaster, we will see that by working in appropriate 
function spaces, and taking advantage of the nature of the nonlinear term, this factor will 
disappear. Its presence is in part due to the fact that we chose to work in "momentum" space, 
rather than "position" space, because the linear problem is most naturally studied in Floquet 



Nonlinear Stability 



12 



variables. If we rewrote these terms in position space (i.e., in the original (x, t) variables), they 
would look much less singular. 

We will work in Sobolev spaces, and we define 

H q , r = {v\(l-d 2 p y/ 2 (l+pY 2 veL 2 }, (3.1) 

equipped with the corresponding norm || • \\ q ^ r . The function w c will be an element of H q r . 

The function u; s has two components. The first component comes from the central branch 
of the spectrum of the linear operator (1.9), and will also be in H q r . The second component 
comes from the stable branches of the spectrum, and it depends on both p, and x. It will be an 
element of the space: 

H q ,r,v = {w = w(p;x) \ w(p;x) = w(p;x + 2n), 

(1 - dlr' 2 {\ - d 2 p Y' 2 (\ +pY 2 W E L 2 (R X [-7T, 7T])} . 



By a slight abuse of notation, we will denote by \\w \\ nq the sum of the H q r norm of the first 
component of w* and the H qrv norm of the second component, and by || w s || q r v , we will mean 
the H qrv norm of just the second component. We will also use H q r v to denote the space of 



all functions with finite H n r „ norm. 



The non-linearities satisfy the following bounds: 
Proposition 3.1. For every q > 2 and every r > there is a constant C for which 

IK/>,e-^ 2 )|| g _ lir < cimi 2 ^ 



,r ) 



||eV 3 (^,e-^ 2 )|| g)r < C||«;|| 3 9)r , 

for all t > 0. 

Proposition 3.2. For every q > 2 and every r > there is a constant C for which 



(3.2) 



|e T / 4 («/> S e- T/2 ,e- r/2 ) 



I q,r 



< Ce^ 2 W\\ Hq ^ {e-^ 2 \\wX, r + e-W\\ Hq J 

x(l + e-T/ 2 \\w%, r + e-T\\ W s\\ Hq J, (3.3) 
\e^ 2 g(w c , w s e-T/ 2 ,e-/ 2 )\\ nq r [/ 

<Ce-(e-^ 2 \\w% !r + e-\\ w ^ Hq J 2 

x(l + e-^ 2 \\w% !r + e-\\w>\\ nq J, (3.4) 



for all t > 0. 

Remark. Note that every factor of 1 1 w c \ \ q jT . is multiplied by e ~ T I 2 and every factor of 1 1 \ \ Hq 
is multiplied by e~ T . 



. r, v 
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Remark. As we pointed out above, the nonlinear terms depend on the constant t . However, 
the bounds in the two preceding propositions are independent of t . More precisely, for any 
T > 0, the the constants C in both propositions can be chosen so that the estimates in (3.2)-(3.4) 
hold for all t > T. 

The proofs will be given in Appendix A. Note that one loses a power of p in the first 
estimate of Eq.(3.2), but of course, one "gains" the square of the function. 

We will regain the "lost" power of p by examining in detail the semi-group generated 
by L. We denote by P N the projection onto the space spanned by the N eigenvalues {u ■ = 
— j /2}j=o,...,N-i °f We define Q N = 1 — P N . (We verify in Appendix B that these 
projections are defined.) On the space corresponding to Q N , we expect the norm of the semi- 
group generated by L to decay like exp(r / u Ar ). This is indeed the case. 

Theorem 3.3. For every e > 0, there is a constant N and a function r(N, q) such that for 
every N > N , every q > 1 and every r > r(N, q), there is a C = C(q, r, N) < oo such that 

V L Q N v\\ qr < (3.5) 

q ' v a ( r ) 

where a(r) = 1 — e~ T and L = —p 2 — \pd p 
The proof will be given in Appendix B. 

We also need an estimate on the linear evolution generated by M exp( -_ T / 2 ). Let U T be the 
solution of 

e~ T d T U T = M exp{ _ T/2) U T , 

with initial condition U Q = 1. (Compare with the linear part of (2.9).) Then, in Appendix B, we 
prove 

Theorem 3.4. If w G H q r v , then there exists c > 0, such that for all r > 0, 

\\U T w \\ ng ^ < exp(-e Cor / 2 )|K|| W9 ^ . 

With the help of the bounds Proposition 3.1-Theorem 3.4, we can now reformulate the 
problem in terms of invariant manifolds. The Eq.(2.9) can be written as an autonomous system 
by defining r] = (t + t )~ i/2 = e" r/2 : 

d T w c = Lw c + V - 2 (f 2 (w c , V ) + f 3 (w c , V ) + f 4 (w c ,w\ V )) , 
r] 2 d T w s = M v w s + r]~ l g(w c , w s r], rj) , (3.6) 
d T V = ~h ■ 

We will construct an invariant manifold tangent at the origin to the eigenspace correspond- 
ing to the largest eigenvalues of L, and the r\ direction. We subdivide the center variable w c 
according to the projection Q N defined earlier, where is fixed once and for all. Define 



x 1 = (1 - Q N )w°, x 2 = Q N w c , x 3 = w s 



(3.7) 
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Note that the variable x l is in a finite dimensional space, while x 2 and x 3 are in infinite 
dimensional Hilbert spaces. The system of equations Eq.(3.6) now takes the form 



d T x i = A x x x + N l (x l ,ri,x 2 ,x 3 ) , 

(3.8) 



d T x 2 — A 2 x 2 + N 2 (x l , rj, x 2 , x 3 ) , 
rj 2 d T x 3 = A 3>r] x 3 + N 3 ( Xl ,rj,x 2 ,x 3 ) . 

Here A l = (1 - Q N )L, A 2 = Q N L, and A 3 r) = M r 

Remark. In view of later developments, we consider x x and rj to be the "interesting" variables 
and x 2 and x 3 the "slaved" variables, hence the new order of the variables. 

Remark. Eq.(3.8) is a very singular perturbation problem, because of the factor of rj in front 
of the derivative of x 3 . What is more, since rj(r) = e~ T ^ 2 , it becomes steadily more singular 
in precisely the limiting regime in which we are interested. Nonetheless, we will see that the 
invariant manifold theorem provides just the tool we need to understand this limit. Singular 
perturbation problems of this type do not seem to have been studied much, but they do arise 
naturally in other contexts, such as the study of parabolic equations in cylindrical domains 
([W2]). 

We shall call Eq.(3.8) the full system. To simplify the notation, we shall omit the depen- 
dence on rj in A 3 . Consider the spectra of A x , A 2 , A 3 . From what we have seen earlier, we 
find that 

spectA) = {0,-l/2,-l,...,-(iV-l)/2}, 
spec(A 2 ) C [-ex), -N/2] , (3.9) 
spec(?7" 2 A 3 ) = [-oo, -c/rj 2 ] , 

where c is some positive constant. Thus, we expect to apply a pseudo center manifold theorem 
to "slave" the variables x 2 , x 3 to the variables x x and rj. While there are certain technical 
difficulties associated with the very singular perturbation, in Appendix C, we demonstrate the 
following Proposition: 

Proposition 3.5. Fix N > 0. There exist r > 0, q > 1, and v > 1/2, such that the system 
of equations (3.8) has an invariant, N + I -dimensional manifold, given in a neighborhood of 
the origin by the graph of a pair of functions 

h* : K N x R - H q ^ r , 
h 3 : ~R N x R -» 1-L q r v . 

We next turn to the task of showing that the invariant manifold we found for Eq.(3.6) 
actually attracts solutions at an exponential rate. 

Notation. It is useful to introduce the notation £ = (x 1 ,rj) for the two relevant variables. 
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Consider a solution of the form (u; c (t), w*(t)) of Eq.(3.6), with w c (t) = (x 1 (t), x 2 (t)) 
as in Eq.(3.7), and w s (t) = x 3 (t). We wish to show that 

(£(r), x 2 ( T ), x 3 (r)) — (£(r) , /^(£(r)) , ^(£(r)) ) , 

as r — > oo, and furthermore, that it does so at an exponential rate, given essentially by the least 
negative eigenvalue, n N , of the operator A 2 . 

Proposition 3.6. Fix N > 0. For every 5 satisfying < 5 there is an e > such that if 
the solution ofEq.(3.6) remains in a neighborhood of the origin of size e one has the following 
bound: There is a C* < oo for which 

\\x 2 (r) - m(r))\\ q ,r + - ^(e(r))||„„ p>1 , < C*e~^-^ , 



as r — > oo. 

Proof. This proof is relatively standard, see e.g., Carr [C]. Let 

z(t) = 

Then we have 



* 2 (t)-/**(£(t))\ /* 2 (t) 



i=| ^ + ^,%^) , (3 . 10) 
where, with the notation of Eq.(3.8), 

n^z^z,) = Nj (tz 2 + h;(0,z 3 + hU0) - ^.(£,^(0,^(0), 

for j = 2, 3. The only novelty in Eq.(3.10) w.r.t. [C] is the factor of r\ in the "3"-component 
which is the reason for our repeating his arguments. But we can integrate Eq.(3.10) explicitly 
and get 



Z 2 



(r) = e rA *z 2 {0) + [ T d*e^-^N 2 (Z(a),z 2 (*),z 3 (<j)) , 
Jo 

— 1/2 

We assume rj(0) > 0, since we are interested in the case 77(0) = t , and we have chosen the 
scaling factor t to be a positive, finite constant. Note also that £ remains in a neighborhood 
of the origin, as r — > oo. From the bounds on the non-linear terms we see that if the solution 
satisfies 

\\x 2 (r)\\ q , r + \\^{r)\\ nqrv < p, 
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for all t > 0, then, with v N = N/2, the modulus of the N th eigenvalue fj, N of L, we have 

\\z 2 (r)\\ q , r < e—\\z 2 (0)\\ qjr + Ce f dae~^>- {\\z 2 {a)\\ q , r + \\z 3 (a)\\ Hq p J , 

Jo 

W^(r)\\n q ^ < e^ r2 -^ 0) -> N \\z 3 (0)\\n q ^ 

(3.11) 

In deriving these inequalities, we used the inequalities 

\\e TA2 N 2 (^z 2 ,z 3 )\\ q , r < e-™»\\N 2 (t,z 2 ,z 3 )\\ q _ hr , 
\\e^N 3 (^z 2l z 3 )\\ ng ^ < e-^\\N 3 ^z 2l z 3 )\\ ng ^, 

which follow from the bounds of Appendix B. If we now fix 5 > and define 

C 2 (r) = sup e->-- 5 )||^ 2 (r')|| 9)r , 

0<t'<t 

C 3 (t) = sup e T '^-^\\z 3 (r')\\ Hq _, 

0<r'<r 

then the Eq.(3.1 1) leads to the inequality 

C 2 (r) < K l+ K 2 e(C 2 (r)+C 3 (r)) f da e~^ s , 

Jo 

C 3 (t) < K 3 +K 4 e(C 2 (t)+C 3 (t)) £ da -^e^r 2 -v(<rr>N e (T-a)(v N -5) 

If we insert into these integrals the definitions 

rj(a) = exp(-a/2)?7(0) , i](r) = exp(-r/2)r](0) , 

we find that both integrals are uniformly bounded in r > if ?y(0) is in a compact subinterval 
of (0, 1). The proof of Proposition 3.6 is complete. 

Thus, all solutions near the invariant manifold approach it exponentially fast in r. 

One can now show without difficulty that every solution approaches exponentially quickly 
a particular solution on the (approximate) invariant manifold 

[x x {t),t] = 0, h* 2 {x x {r), 0), /^(^(r^O)) . 

This consists simply in translating the pp.2 1-24 of [C] into the present setting and thus there is 
no need to repeat this argument here. 

If we combine these results with Proposition 3.5, we arrive finally at a description of the 
invariant manifolds which exist close to the origin for (3.8). 
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Theorem 3.7. Fix N > and 5 > 0. There exist r > 0, q > 1, and z/ > 1/2, such 
thai: the system of equations (3.8) has an invariant, N + I -dimensional manifold, given in a 
neighborhood of the origin by the graph of a pair of functions h\ : R^ x R — > H q r , and 
h 3 : TL N x R -> H q r il . Any solution of (3.8) which remains in a neighborhood of the origin 
for all t > approaches a solution of the N + 1 -dimensional system of ordinary differential 
equations 

d T x 1 = A lXl + N 1 (x 1 ,rj, h^x^rj), h^x^rj)) , 
<9 T r/ = -±77 , 

which results from restricting (3. 8) to this invariant manifold. Furthermore, the rate of approach 
to this manifold is C(exp(-r(iV/2 - 5))). 

Remark. This theorem almost suffices to prove Stability Theorem 1.2. In particular, it 
emphasizes that in a neighborhood of the periodic solutions of (1.1) there exists a family of 
invariant manifolds, M 2 , M 3 , . . ., described in that theorem. The one remaining piece of the 
puzzle is to describe the behavior of solutions restricted to the invariant manifold, and that we 
do in the next section. 



4. The projection of the non-linearity onto zero momentum 

We have already shown that there exists a (smooth) invariant manifold, parameterized by 
(£ 7 ^(Oj ^(O)' where £ = (x 1 , rf). This manifold satisfies the equation Eq.(3.8), which, in 
the case of N = 1, i.e., in the case of a two-dimensional invariant manifold amounts to 

d T x x = ^(^,77,^(0,^(0) , 

^="^' (4.1) 
<9 T (/i*(0) = A 2 hm + N 2 ( Xl , V ,h*(0,h^)) , 

V 2 d T (m)) = A 3 h;(0 + N 3 (x vVl h*(0,h* 3 (0) ■ 

Note that because iV = 1 the operator A { equals zero (which is the highest eigenvalue of L). 

To understand the dynamics inside this invariant manifold, we now state and prove the 
following proposition, which is based on Schneider's beautiful observation: Let N l (x ll i]) be 
the r.h.s. of the first equation in (4.1), i.e., d T x x = N Y (x Y ,ri). 

Proposition 4.1. There is an x l > such that N l (x l ,0) = 0, for all \x t \ < x l . 

Thus, the non-linearity vanishes identically at "infinite time," which corresponds to 77 = 0. 
Before proving Proposition 4.1, we show that it implies the following important 

Theorem 4.2. If x l (0) is sufficiently close to 0, then there are a constant C < 00 and an x\ 
such that 

\x x {t) -x\\ < Ce~ T/2 . (4.2) 

Proof. Using the fact that r/(r) = e~ T ^ 2 , we can rewrite the equation for x x as 

d T x x = N^e^' 2 ). (4.3) 
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Since N l is a smooth (at least C l+a ) function with N l (x l , 0) = in some neighborhood of the 
origin, there exists a constant C N > 0, such that \N 1 (x 1 , e~ T ^ 2 ) \ < C N exp(—r/2), for 
sufficiently small. Integrating (4.3) and applying this estimate yields: 



\ Xl (r f ) -x^rJI = 



f rfaiV 1 (x(a),e-' T/2 ; 



< dae~ a/2 = 2C Jv e- Tl/2 (l -e (Tl - Tf)/2 ) 



This estimate immediately implies the behavior claimed in Theorem 4.2. 

Proof of Proposition 4.1. The basic idea is to relate iV^Xj, 0) to the non-linear term of 
another problem, which is known to be 0. This other problem is the center manifold equation 
for the perturbations of a stationary solution of Eq.(l.l) restricted to a space of lit -periodic 
functions. In this case, the equation analogous to Eq.(1.12) is 

where F(v) collects the non-linear terms in v. The spectrum of L per is pure point, with a 
simple zero eigenvalue, and all others negative, and bounded away from 0. The eigenvector 
with eigenvalue is u' e , where u e is given by Eq.(1.5). If we call x l per the coordinate in the 
u' e direction, then there exists a one-dimensional center manifold, tangent to this direction and 
given as the graph of a function H(x 1 per ). A very nice observation by Schneider is that this 
center manifold must coincide with the translates of the stationary state u e , which is formed 
of fixed points of the Swift-Hohenberg Eq.(l.l). Hence, on this center manifold we must have 
x l pei = 0. Using this information, the equations for this center manifold take a particularly 
simple form. Let P per denote the projection onto u' s and let <5 per = 1 — -P per . Then the preceding 
discussion implies that the flow ip t per is the identity on x l per , and hence the equations for the 
invariant manifold read: 

*l,per = ^Kper^Kper)) = 0, (4.4) 
,0 

H(x hpsi ) = / rfre-^ L -Q per F(x 1;per ,if(x liper )) 

J — oo 

= -(Qper^per)" 1 ^^^!^, ^(^l,per)) • (4-5) 

We now wish to use this information to prove Proposition 4. 1 . The rough idea is to show that 

N^x^O) = P per F(x liper ,^(x liper )) , (4.6) 
and this quantity vanishes by Eq.(4.4). More precisely, we shall show: 

Proposition 4.3. The cubic term in x i of N l (x 1 ,r]) coincides in the limit 77 — > with the 
cubic term in x l of P F(x l per , H(x l pei )). All other terms in N l go to as 77 — > 0. 

Remark. Since P peT F(x l peT , H(x l pes .)) = 0, this proves Eq.(4.6) and thus Proposition 4.1. 
Proof. The proof of Proposition 4.3 will be given in Appendix D. 
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5. Completion of the proof of Stability Theorem 1.2 

We now consider exactly how the results of the previous two sections about the behavior of 
solutions in, and near, the invariant manifold translate back into statements about solutions in 
terms of the original variables. We will focus specifically on the case considered in the previous 
section in which the invariant manifold is two-dimensional, with coordinates (x 1 ,rj), but the 
results can be immediately extended to the case of a manifold of arbitrary dimension. 

Suppose we have a solution w T = w c T + w s T , of the system (3.6), which remains in a 
neighborhood of the origin for all r > 0. This will be the case if its initial condition is 
sufficiently small in H r © H q r „. We measure the size of w in the norm ||| • |||, which is the 
sum of the H n _ norm of w c , and the H n _ „ norm of w s . By the results of Theorem 3.7, we 
know that there exists a solution, w™, on the invariant manifold such that 



\w r -w™\ 



< Ce- T{l ' 2 - S \ (5.1) 



with 5 > 0. In addition, from Theorem 4.2, we know that there exists some w*, which lies in 
the invariant manifold for which 



\w 



™-w*\\\ < Ce~ T/2 . (5.2) 



Here, w* is the function whose coordinates in the invariant manifold representation is just the 
limiting point x * in Theorem 4.2, i.e., w* = (x\ , 0, h^ixl , 0), h^(x* , 0)). Combining (5.1) and 
(5.2), we see that for solutions that remain near the origin, there exists a function w*, for which 



— w 



< Ce-^l 1 -^ . (5.3) 



Our final task is now to untangle the various changes of variables which we made in the 
original equation. If we first "undo" the rescaling in (2.8), we see that the solution v(£,t), 
corresponding to w(-, r) = w T is 



v(£,t) = w c (sign(£)^\A £ \(t + t ),log(t + t )) 



+ (i + to y /2 w s ^gn(£)^\A e \(t + t ), log(t + t )) (5-4) 
= v c (£,t)+v\£,t). 

One can make a corresponding decomposition of v*, the solution corresponding to w* . 
First consider v c . From (5.3), one has 

K-<' c n; iP = J d P \(i-d 2 p y/ 2 (i+pY 2 (w c (p,T)-w*> c (p,T))\ 2 < ce-^~ 2S \ 

(5.5) 
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According to (5.4), w c (£, r) = v c ($~ 1 (pe~ T / 2 ),t), so substituting this expression — and the 
analog for w* ,c — into (5.5) one finds that the left hand side of that inequality is equal to: 

J dp\(\-dlr/\\+py'\v c ^-\pe-^ 2 ),t)-v^- l (pe-^ 2 )^))\ 2 

> J dp\{\ + p 2 ) q / 2 {v c {^-\pe- T / 2 ),t) - v^ c {^-\pe- T/2 ),t))\ 2 (5.6) 

where in the last integral we changed the integration variable to £ = <E> _1 (pe~ T ^ 2 ) = ^~ 1 (p(t + 

g- 1/2 )- 

Remark. We dropped the derivatives with respect to p in the second line of (5.6) for 
simplicity — one could retain them at the expense of complicating the following expressions. 

Since $(#) x, for x small, and is equal to a constant times x for \x\ large (due to the 
definition of A e ), we see that combining (5.5) and (5.6) and recalling that t > 0, one finds: 

J Ml(l+m ^ () _ v , W < C( -3«,- M ). ( , 7) 

Analogous estimates hold for the "stable" part of the solution. Proceeding as above, one 
can show that 

J](l + 72 2 r / d£\(l + £ 2 ) q/2 (v s (£,t) - v*' s (£,t))\ 2 < Ct~ 5/2{l - 25) . (5.8) 

Thus, the "stable" part of a solution near the origin approaches the solution v* on the invariant 
manifold faster than the "center" part of the solution. (An effect that is entirely in accord with 
one's intuition.) 

We next take a closer look at the solution w* (or v*) on the invariant manifold. From the 
computation in the previous section, we know that since the eigenfunction in the x l direction is 
exp(— p 2 ), cf. Eq(2.14), we have w*(p) = c* exp(— p 2 ) + h^(c* exp(— p 2 )). If we now rewrite 
this in terms of the v(£, t) variables, we find 

v*(i,t) = c*e~ Ait + t- [ / 2 h;(c*e- Ait ) . (5.9) 

Thus, if v(£, t) is a solution of (1.12) (in the unsealed variables), we see from (5.7)-(5.9) that in 
the L 2 ((l+£ 2 ) q / 2 d£) norm, 

v(£, t) = c*e~ Alt + 0{r l/2{l - 2&) ) . (5.10) 

But we know from Section 2 that A e = £ 2 + 0(£ 3 ) for £ small, and A e = c£ 2 , for \£\ large, so 
one finds by an easy and explicit estimate that 

J d£\(l+£ 2 ) q / 2 (e- A * t -e- e2t )\ 2 < Cr 1 ' 2 . (5.11) 
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Combining (5.10) and (5.1 1) one has 

Proposition 5.1. If v is a solution of (1.12) with sufficiently small initial condition (in 

H q,r © H q,r,v)> then 

( J d£\(l + fy/ 2 (v(£,t) - c*e- ft )\ 2 ) 1 / 2 < Ci" 1 ^ 1 " 25 ) . 

Note that if we transform back to the (x, t) variables, this implies the asymptotic estimate 
in Stability Theorem 1.2, and hence the proof of that theorem is complete. 

A. Bounds on the non-linearities 

In this section, we prove Proposition 3.1 and Proposition 3.2. We begin by studying the kernels 
K 2 (£, k), and K 3 (£, k) introduced in (2.7). 

Lemma A.l. There is a constant C such that 

\K 2 (£,k)\ < Cemax((\k\ 2 + \£\ 2 ),l) . 

Proof. By the definition of Eq.(2.7), we have 

K 2 (£,k) = J dxTp e (x)u £ (x)cp k (x)cp £ _ k (x) . (A.l) 

Since u £ and cp k are both uniformly bounded, we have immediately that \K 2 (k, £) \ < Ce. The 
crucial observation of Schneider[Sch] is that because of Eq.( 1 . 1 1 ), repeated here for convenience 

<PeA X ) = U 'e^ + U 9e( X ) + ' ( A " 2 ) 

(with real g e ), K 2 has an expansion 

J dxu e {x){u' e {x)f + u £ (x)(u' £ (x)) 2 (-i£ + ik + i(£-k))+eO(£ 2 + k 2 ) . (A.3) 

Note that the first term vanishes because u is a symmetric function and hence u(u') 3 is odd, and 
the term which is linear in k and £ vanishes as well, because of momentum conservation, so the 
proof of Lemma A. 1 is complete. 

Remark. Note that a similar calculation immediately shows that the kernel K 3 satisfies: 

1^3(^1^2)1 < Ce. 
We now need the following auxiliary result: 
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Lemma A. 2. If p 2 and p 3 are in H q r , and if p l = p l (p, p) is a C r function, then 

H(p) = J dp p 1 (p,p')p 2 (A(p,p,r))p 3 (p') 

is in H n „ and 

IMI^r < C||Plllc^llP2llq,rllP3llq,r • 

Proof. Recall from Eq.(2.12) that A(p, p', r) ss p — p , so we are really estimating a slightly 
distorted convolution. If A(p,p', r) = (p — p'), the proof is easy using the definition of the 
norms. In the present case, where A(p, p', t) is not trivial, the result follows in a similar way by 
"undoing" part of the variable transformation which led from the variables £, k to the variables 
p, p . To simplify matters, we consider only the somewhat easier problem of bounding 

J dp' <S>'(pe-T/ 2 )p 2 (A(p,p',T))p 3 (p') . (A.4) 

Using the definition of A(p, p', r) this is equal to 

J dp'^{pe- r / 2 )p 2 (e T / 2 ^- 1 (^{pe- T / 2 ) - ^^e"^ 2 ))) p 3 (e^" 1 ^(p'e"^ 2 ))) . 

(A.5) 

Changing variables to k = <E»(e T ^ 2 p) and £ = <£>(e T ^ 2 p'), we get 

J dle T/2 p 2 (e T/2 $-\k - £)) p 3 (e T/2 <f>- l (£)) . (A.6) 
We now define a function \T/ T by 

# r (e r/ V) = e T/2 ^~\x) , 

and note that from = x • (1 + O(x)) it follows that * T (y) = y • (l + C(e" r/2 y)). We 
can rewrite Eq.(A.6) as 

J d£e^ 2 p 2 ^ T {e^ 2 {k-£)))p 3 ^> T {e^ 2 £)) . (A.7) 

We define next pAk) = p- o \l/ r , and we see that Eq.(A.7) is equal to 

J d£ p 2 (k - £) p 3 (£) . (A.8) 

Thus, we can bound the H r norm of Eq.(A.4) by ||p 2 || g r ||P3ll g r , and, since \T/ T is uniformly 
close to the identity for all r, this is in turn bounded by const. ||P2llq,rllP3llq,r- This proves 
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Lemma A.2 in this special case. The extension to the general case is easy and is left to the 
reader. 

We now have the necessary tools to attack the proofs of Proposition 3. 1 and Proposition 3.2. 

Proof of Proposition 3.1. If we write out the transformation leading to f 2 , i.e., from 
Eq.(2.7) to Eq.(2.9), we get, using Eq.(2. 10), 

e T • {f 2 (w,e-^ 2 ))(p) = e-3x($(pe-/ 2 )) f ' dp' e~^' (p' e~^ 2 ) 

J-P(r) (A.9) 

xK 2 ^{pe- T / 2 )Mp'e- T/2 ))w{A{p,p\r))w{p'), 

where 

P(t) = *-\\)e T ! 2 « \e Tl2 . 

Webound|K 2 ($(pe- r/2 ),$(p'e- r/2 ))|by Ce|$(pe _T/2 ) 2 + $(p / e _T/2 ) 2 | using Lemma A.l. 
Since the expressions $(pe~ r / 2 ), and <&(p'e -T ' 2 ), in Eq.(A.9) are bounded, and <&(x) = 
x(l + O(x)), we can extract another factor of e~ T / 2 and get a bound on e T / 2 of the form 

const. e r/2 x($(pe- r/2 )) / dp' (|$(p e - r/2 )| + |$(p' e - r/2 ) |) • \w(A(p,p' ,r))w(p')\ 

J-P(t) 

/oo 
dp' \p + p'\\w(A(p,p', t))w(p')\ . 
-oo 

(A. 10) 

If w is in H , then with the aid of Lemma A.2, we can estimate the H q _ l r norm of Eq.(A.9) 
by C||mj|| 2 r . Note further, that from the above discussion it is also clear that e T f 2 (w c , e~ T ^ 2 )(p) 
is also a smooth function of e~ T ^ 2 . 

Remark. The factors \p\, \p'\ are responsible for the loss of one power in the norm estimate 
of Proposition 3.1. It is only in the study of the flow within the invariant manifold that we will 
need the second order bound of Lemma A.l. 

Remark. Note that the nonlinear terms depend (implicitly) on the constant t which entered 
the definition of the new temporal variable r. However, all the estimates above (as well as those 
which follow in the proof of Proposition 3.2) are independent of this constant. 

The bound on / 3 is similar, but no additional regularization is needed, since there are two 
integrations, each of which contributes a factor e~ T ^ 2 . We leave this to the reader. The proof of 
the asserted bounds of Eq.(3.2) is complete. 

We now turn to the estimates of the nonlinear terms / 4 and g. Because these terms involve 
the w s , we begin with a discussion of the appropriate function space for these components. These 
were defined in Section 3, but we repeat them here for convenience. Recall that w c G H q r , 
while w s G H q r © H q r u , where 

Hq,r,u = {w = w(p; x) \ w{p; x) = w{p; x + 2n) , 

(1 - d 2 r /2 (l - a 2 ) r/2 (l + P 2 ) q/2 W G L 2 (R X [-7T, 7T])} . 
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The fact that w s is an element of the direct sum of two spaces reflects the fact (see the paragraph 
preceding (2.6), and then (2.8) ) that it has two components, the first of which comes from 
the central branch of the spectrum of L e , but with £ localized away from zero, and the second 
component coming from the stable branches of the spectrum of L £ . In a slight abuse of notation 
we will denote by ||w s || w the sum of the H q r norm of the first component of w* and the 
H qrv norm of the second component, and by ||w s || g)T . ^ we will mean the H norm of just 
the second component. 

Remark. An easy fact which will be useful later is that if we expand w(p; x) G H in a 
Fourier series with respect to x, 

oo 

w(p;x) = einX <(p) > 

n= — oo 

then the H qru norm of w is equivalent to the norm 

oo 

iMiiw = E a+^riKiilr- (a.ii) 

n= — oo 

Thus we will use the two norms interchangeably. 
Now consider 

e T f 4 (w c ,w s e- T/2 ,e- T/2 ) . (A.12) 

We shall concentrate on the most "dangerous" piece which is the quadratic term with one factor 
of w c and one of w s . Other terms are "less dangerous" in the sense that they contain either more 
factors of w* each of which contributes a small factor of e _T / 2 , or more convolutions which 
again contribute a factor of e~ T I 1 . The quadratic piece of (A.12) has the form 

e r 3x($(pe" r/2 )) J dx<p^ e - T/2) (x) u £ (x) 

/ p ( T ) 
dp' e" r/2 $ Ve~ r/2 K (A(p, P, r)) 13 ^ 
-P(t) 

x < Pr(p,p>,T)( x ) e ~ T/2w *(p'> x ) ■ 

As we mentioned above, w s has two components — one in H , and one in H . The 
contribution from the component in H is bounded by the same techniques used to control 

/ 3 — note that it is not necessary to extract any additional factors of e~ r / 2 , since we get one from 
the integration, and one from the fact that each factor of w s is multiplied by e~ T//2 . Thus, we 
restrict our attention to the component of w s in H r v , which is where the new ingredients are 
necessary. 

Interchanging the order of the x and p integrals, we use Lemma A. 2, with 

Pl (p,p') = sup|3x($(pe- T / 2 ))$ / (p'e- r / 2 )^ (pe _ T/2) (x)^ r(pp , )T) (x)w e (x)| , 

X 

p 2 (r) = \w c (r) \ , 
Piip') = | / dxw s (p';x) \ . 
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Since <Pt(x) and u £ (x) are smooth, 27r-periodic functions of x, and ||pi||cr is bounded, the 
Lemma A. 2 implies that the H norm of (A. 13) is bounded by 

C|KIUI J dxw\--x)\\ q ^ r . (A.14) 

The H q r norm of the integral can be bounded by 

supKMIU < C\\w*\\ Hqr ^ (A.15) 

X 

provided v > 1/2, where we used Sobolev's inequality to estimate the supremum over x. 
Inserting (A.15) into (A.14) yields the bound claimed in (3.3). 

The remaining terms in / 4 can be bounded in a similar fashion, but as noted above, they 
will tend to as r — > oo. In fact, they will be bounded by Cee~ T ^ 2 . 

Proof of Eq.(3.4) of Proposition 3.2. We finally bound the non-linear term 

e T/2 g(w c , w s e~ T/2 , e" r/2 ) . (A.16) 

In bounding e T ^ 2 g(w c ,w s e~ T ^ 2 ; e~ T ^ 2 ), recall that just as did, this expression will 
have two components — one in H , and one in H . The component in H is bounded 
using exactly the same techniques used to control the term f 4 above, so we concentrate here on 
explaining the new ingredients necessary to bound the component mH qru . 

As in the bound on / 4 , the potentially largest terms are those of minimal order, because 
each additional order provides a factor of e~ T ^ 2 . So we look at the terms which are quadratic 
and which are of order w c w c , w c w s , and w*w s , respectively. The first term leads us to study 

e T ' 2 P^(u £ {x) ^^V c (p-p^ C (P>p'(^)^-p'W) • (A.17) 
Rescaling as in (2.8), we see we must bound 

Pi {pe -r n {u £ (x) dp'&ip'e-^ 2 ) 



x 



P*(p'e— /^(s) Vripy ,t)( x ) w °(p) w c (A(p,p\ r)) 



(A.18) 



Note that the prefactor of e T ^ 2 has disappeared due to the factor of e~ T I 2 which we gain as usual 
from the change of variables. 

Since the projection P e has bounded norm and is a smooth function of t, we can discard 
this factor at the price of introducing an overall constant in the estimate. Note next that the 
square of the H q r u norm of the remaining expression is equal to: 

• P(t) 

dp' &{p'e- T/2 )w c {p')w c {&(p,p', t)) 
" yJ W (A. 19) 



X 



\ u e( x ) ( P®(p'e-r/2)( x ) ¥t{ p ,p',t)( x ) llff(dx) 



2 

H q , r (dp) 
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where the H v norm is the if^-Sobolev norm of the quantity 

U e( X ) ¥>*(p'e— /2)(a0 Vr(p,p',r)(^) , 

considered as a function of x, and the H r norm is the norm of the resulting function of p. Since 

u £ {x) <p<f>(p> e -T/2)(x) ^r( PP ' t)( x ) * s a smooth function of x, p' , and p, there exists a smooth, 
bounded function ip(p, p'), such that 



(A.20) 



V>(p,p') = ll« e (x)^ (p/e _ T/2) (a:)v7r( P , P ',r)(^) lln^ete) • 
But now, by Lemma A.2, we can conclude that (A. 19) is bounded by 

|| / P(T) ^$'(pV-/ 2 )«; c (p> c (A(p,^r))^(p,p')l| 2 H 9r ( dp ) < C||*Vllc-Klll,r- 

J-P(t) 

(A.21) 

We next consider the quadratic term in g which contains one factor of w c and one factor 
of w s . In this case, the analog of (A. 18) is 



-Hp 



u e (x) 



P(r) 



dp'&(p'e- T/2 ) 



x 



-P(t) 

^(p.p'.r)^) w c (A(p,p',T))w a (p' ;x) 



(A.22) 



Note that in this case, we pick up an extra factor of e~ r ^ 2 , in comparison with (A. 18), since 
each factor of w* is multiplied by this exponential. 

Once again, we must contend with the fact that w s has two components. However, the 
component in H r behaves exactly as in the estimates leading to (3.2), so we concentrate on 
the component mH qrv . 

As above, the projection operator can be dropped at the cost of an overall constant, and we 
are left with the task of bounding the H q r u norm of the remainder. The square of this norm is 
equal to 



rP(r) 

/ dp'<f>'(p'e- T/2 )w c (A(p,p',T)) 



x \\u e (x) ¥>*(p/e-T/2)(aO (p r(pjptjT) (x) w s (p'; x)\\ Hv{dx) 

< c\W\\ cr \\w c f Har 



H q , r (dp) 



(A.23) 



X 



K( X ) ^(p'e-T/i) ( X ) <Pr(p,p>,T) ( X ) W °(P > X ) llff(dx) 



H q , r (dp') 



by Lemma A.2. Note that the pair of norms on the last factor is equivalent to computing the 
square of the H qru norm of 



u e( x ) <P*(p> e -T/*)( x ) ¥r( P ,p',T)( x ) w s (p';x) . 



(A.24) 
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Since u e , <£>, ip r , and A are all smooth, bounded functions, we see just by writing out the 
definition of the norm that this is bounded by 

C\\w s \\ 2 Hq ^ . (A.25) 

If we estimate the term quadratic in w* in a similar fashion, and combine this estimate with that 
in (A. 21) we see that the quadratic terms in e~ T ^ 2 g(w c , w s e~ T ^ 2 ; e~ T ^ 2 ) are bounded in H q r u , 
by 

C{\\w c \\ Hqr +e-^ 2 W\\ Hqr J 2 . (A.26) 
Analogous estimates of the cubic terms lead to a bound 

Ce-^dl^H^+e-^Hu;'!!^) 3 , (A.27) 

where the additional factor of e~ T l 2 comes from the additional convolution. Combining (A.26) 
and (A.27) leads to the estimate in (3.4) and completes the proof of Proposition 3.2. 



B. Bounds on the linear operators 

In this Appendix, we give bounds on the semi-group generated by L and on the linear evolution 
defined by M exp( _ r/2) . 



B.l. Bound on the semi-group generated by L 

We consider the semi-group whose generator is L = d 2 x + \xd x + \. Note that in this section, 
for ease of use, we define L in the Fourier transformed variables, compared to Section 2. Fourier 
transformation is an isomorphism from H r (in the p-variables) to H r (in the x-variables), so 

establishing estimates on the semigroup associated to d 2 + \xd x + \ in the space H r q (dx) will 
immediately imply estimates on the representation of L in the p- variables in the space H r (dp) . 
In order to avoid confusion, in what follows we will denote by I • L _ the norm on H r Jdx). 
With this notation, the norms || • || and | ■ | resp. the spaces H q r (dp) and H r q (dx) are 
equivalent. 

The integral kernel of the semigroup generated by L is given by [GJ] 

(e TL v)(x) = ; 1 f dze- z2 ^ 4a(T)) v(e T / 2 (x + z)) , 
V47ra(r) J 

where a(r) = 1 — e~ T . If we denote by T the operator of multiplication by exp(x 2 /8) and 
by H Q the harmonic oscillator Hamiltonian H Q = d 2 x — x 2 / 16 + 1/4, (note the unconventional 
sign!), then 

L = T~ l H T . 

Thus, the two operators L and H are "the same," but they act on two quite different spaces. If 
the {<Pj}j >0 are the eigenfunctions of H , then the ^ ■ = T~ x ip- are the eigenfunctions of L, 
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with the same eigenvalues [i- = -j/2. We let P n f = J2j< n V'j(V'j) f) q > where (•, -) q is the 
scalar product 

(f,g) q = (Tf,T(l-L o rg) = (Tf, (1 — H ) q Tg) . 
We next show that for n < q — 2, the operator P n is bounded in H r q (dx). First of 
all, the eigenfunctions (p ■ are bounded by 0(l)\x\^e~ x ^ at large x. Therefore, we also have 
ipj = T~ l ip- e H rq (dx), since it decays exponentially. Finally, 

(^J) q = (T^,(l-H o yTf) = \l-^( Vji Tf) t 

and the last scalar product is bounded if / e H r q {dx) when r > j + 2, since, with a weight 
function W{x) = (l+x 2 ) l/2 , 

K^.,t/)| < c|(^',/)| < ck^- 1 ,^ 1 /)! 

< C||^ +1 /Il 2 < C|/| 0>r . 
Thus P n is defined. We let Q n = 1 - P n (in H r q (dx)). 

Theorem B.l. For every e > 0, there are an m and a function r(m, q) such that for every 
m > m , every q > 1 and every r > r(m, q), there is a C = C(q, r, m) < oo such that 

\e TL Q m v\ qtr <^^e^^\v\ q . hr . (B.l) 

v a ( r ) 

Remark. The function r(m, <?) is of order C(m + (/). 

Proof. To explain the strategy of the proof, we need some notation. Let P^ denote the 
projection in H q {dx) onto the subspace spanned by and let = 1 — P^ . Then, 

formally, TQ n = Q^T, and LQ n = T^HqQ^T. This suggests that L restricted to Q n has 
no spectrum in the half -plane {z | Re z > —\/J, n+i | }, and thus one can understand the decay in 
Eq.(B.l). The square-root singularity at r = is related to our gain in smoothness. The problem 
is that TQ n = Q^'T is ill-defined. However, it will be well defined if we localize near x = 0. 
In that region, the heuristic argument will be seen to be valid, whereas in the complement of 
such a region, when \x\ > R, decay will be shown by direct methods, using the explicit form of 
the integral kernel. 

We study first the quantity X fl e rL , where Xr is a smooth characteristic function which 
vanishes for \x\ < R and is equal to 1 for \x\ > 4R/3. Thus we study a region far from the 
origin. Our bound is 

Proposition B.2. For every q > 1 and every r > there exists aC(q,r) < oo such that for 
all v E H r q (dx) one has 

\XR^ L v\ q!r < ^ e -/ 2 (e-/ 2 + e -3tf , (B . 2) 
y a \ T ) 

\x R e TL v\ q , r < C(q, r)e^ 2 ( e — / 2 + e " 3rf / 16 ) ^ . (B3) 
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Corollary B.3. For every q > 1 and every r > there exists aC(q,r) < oo such that for all 

v G H rq (dx) one has 



V ft ( T ) 

le^l^ < C( ? ,r)e^/>| g;r . (B.5) 



Remarks. The improvement over [W] is that we "gain" a derivative in x. The corollary 
follows easily by repeating the proof of Proposition B.2 with R = 0. 

Proof. We let D = d x and denote, as before, by W the operator of multiplication by ( 1 + x 2 ) 1 ^ . 
Then 

\x R e TL w\ 2 qir and ^ U ' ; ' - 

<?'<<? 



are equivalent. We shall only consider the term with the highest derivative, because only there 
is the issue of regularization important. Thus we are led to bound 



X 2 = \\W r D« XR e TL w\\ 2 2 . 
Since L = d 2 x + \xd x + \, a quick calculation shows that 

We Th = e Tq/2 e rL D q . 
The diverging factor exp(rq/2) will appear in the final bound. Note now that 

(e rL DV) (x) = - 1 e T/2 / dze- z2/(4a{T)) (D q v)(e T/2 (x + z)) , (B.6) 
A/47ra(r) J 

which upon integrating by parts becomes 

Use now the Schwarz inequality in the form (for positive / and g), 

\\f*g\\l = J dx J dz i dz 2 f{ z \)f{ z i)9{x - z x )g{x - z 2 ) 

< J dz l dz 2 f{z l )f{z 2 )\\g{- - z x )\\ 2 \\g{- - z 2 )\\ 2 

\ 2 

dzf(z)\\g(--z)\\ 2 
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This leads to a bound 

,rq/2 



X < -=== / dz-LL e-/( 4 ^))||^ Xi? ( Dg -i w ) (e ^/2 ( . + 2)) || 2 



(B.7) 



^1 + ^2 5 



where we let R Y = {x : \x\ < lR/%} and R 2 = R \ R x . To be more precise, we define by 
the scaling of a fixed function: Xr( x ) = x( x /R)- If -R — > oo, then c^x^x) = C(-R _1 ) and 
therefore it is uniformly bounded. 

Lemma B.4. (Lemma A.2 of [W]). One has the bounds 

Proof of Lemma B.4. Consider first the case \z\ < 7i?/8. Since |x| > i? on the support of 
Xr, we have \x + z\ > \x\/S and hence 

(1 + x 2 )/(l + (e T/2 \x + z\f) < const. e" r . 

Using this, we bound 

/ dx(\+x 2 Y\x R (x)v(e^ 2 (x + z))\ 2 

= \ dx ■ (1 + (e^\x + z\) 2 Y\v(e^ 2 (x + z))\ 2 

Jri (1 + (e T / z \x + z\Y) 

< const. e~ Tr e~ T ^ 2 \v\l !r < const. e~ Tr \v\l r . 

In the second case, we get 

dx{\ + x 2 ) r \x R (x)v(e T/2 (x + z))\ 2 



R 2 



(l + y 2 Y 

< const. e- T ' 2 {\ + z 2 Y\v\ 2 0r < const. (1 + z 2 ) r \v\l r . 

The proof of Lemma B.4 is complete. 

Continuing the proof of Proposition B.2, we first bound the integral over R { in Eq.(B.7). 
We get from the first alternative of Lemma B.4, 



A/47ra(r) Jr x ^2a(r) 

< COnSt . 1 P Tq ' 2 f dz ^_ e - Z V(4a(r)) e -W2 



V 47ra ( r ) ^ V /2a ( r ) 



9-1, r 



1 



< const. - 1 e ^(g/a-r/2) 



9-1, t- • 
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Similarly, using the second alternative in Eq.(B.8), we get 

X 2 = J^r F f q/2 f R ^- 7 ^e- 22 /(4«M)||^ Xi? ( D ,-i w ) (e ^ ( . + 2)) || 2 



< const. - T ^=e rq/1 [ dz(l + z 2 ) r/2 )^ e jv 
A/47ra(r) Jr 2 \j2a{r) 



z 2 /(4a(r))|.., i 

q-\,r 



< COnSt. 1 e r 9 /2 e -3flVl6^, ; 

V 4 ™( T ) 



q-l,r ' 



since 3/16 < (7/8) 2 /4. Note that the constants above depend on r and q, but can be chosen 
uniformly for all i? > 1. The proof of Eq.(B.2) is complete. Omitting the integration by parts 
in Eq.(B.6), the assertion Eq.(B.3) follows in the same way. The proof of Proposition B.2 is 
complete. 

We next study e rL Q n (l — Xr) w - We have the following bound 

Proposition B.5. For every e > 0, q > 1, and every r > there is a C(e, q, r) < oo such 
that 

\e TL Q n (l ~ X R )M q ,r < ^4e-^^e R2 ^\ W \ q - hr • (B.9) 

V a \ T > 

Proof. Recall that T = e x2/8 and that L = T~ l H T. The operator T( 1 - x fl ) is bounded and 
11^(1 - x fl ) || < const. e Rl/6 . Therefore we have 

Q n T{\ - Xr ) = (1 - P n )T{\ - x R ) = T{\ - Xr ) - TP£\\ - Xr ) 
= T(1-PW)(1- Xb ) =TQi°)(l-x B ), 

where is the orthogonal projection onto the complement of the subspace spanned by the 
first n eigenvalues of H in H . It is easy to see that on H (dx), the operator (1 +x 2 ) 1 ^ 2 (l — 
if ) -1 / 2 is bounded. Thus, we get, using the spectral properties of H (on Q^), 

\e TH °TQ n (l- XR )w\ qjr = r" 1 / 2 

x |(1 - H )"V(t(1 - H )y/ 2 Q^T(l - Xr W 
< const. r- l / 2 \e^(r(l - H Q )) l / 2 Q^T{\ - X r)wL_ 1 



q,r 
q-l,r 



< const. r - 1 /2 e - r|M "+ l|+Te |T(l - X R )w\ q -\, 

< const. r- 1/2 e- r| ' i "+ l|+r£ e i?2/6 |«;| 



lg-l,r • 

(B.10) 

The proof of Proposition B.5 is complete. 

End of proof of Theorem B.l. We first rewrite e rL Q n as 

e TL Q n = e^Q^' 2 = e- L / 2 Q nXR e TL/2 + e- L / 2 Q n (l-X R )e TL/2 . 
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The second term can be bounded by Proposition B.5 and Eq.(B.5) as 

V a ( T ) 

< - c R2/6 ~ T ^ n+[l/4 c Tq/4 \w\ 

This quantity is bounded by 

C „-rn/8u„ 

l<j-l,r 5 



g-l,r • 



( e-™/VL-i. P . (B- 11 ) 



provided n is much larger than q and i? 2 / 6 < rn/ 16. The first term can be bounded by Eq.(B.5) 
and Eq.(B.2) as 



V a ( T ) 



q-i,r 



< ^= eT9/2 ( e - Tr/2 + e " 3i?2/16 ) l-l.-i. (B-12) 
C 



-ttz/8|„„| 

<2-l,r ) 



< =e-™/ 8 |«, 



provided r > n/4 + q and 3i? 2 /16 > r(n/8 + <z/2). Note that the conditions on R from the 
first and second term are compatible 

Combining Eqs.(B.ll)-(B.12), we get 

\e TL QM q ,r < -^=e-™/*\w\ q _ hr . (B.13) 

v a ( r ) 

It remains to improve the decay rate from nj 8 to |// m+1 |. The idea is to just take n = 8(m+l). 
Then we find 

e TL Q m = e TL Q n Q m + e Th P m Q m + e rL (P n - P m )Q m . (B.14) 

The first term is bounded by Eq.(B.13), and m/8 > — |/x n+1 |. The second term vanishes and 
the third is diagonalized explicitly: 

e TL (P n ~ PjQm = T- l e-^T(P n ~ P m )Q m = T~ x e^{P^ - P£ ) )TQ m . 

We are operating here on the finite dimensional subspace spanned by the eigenvectors (p m+l , 
. . . , (p n , and there the technique of Eq.(B. 10) yields a bound 

C 



Tl/7 \p-\^rn+l\T 



Combining this with the bound on the first term in Eq.(B.14), we complete the proof of Theo- 
rem B.l. 
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B.2. The linear evolution generated by M v 2 

In this section, we deal with the problem of giving bounds on the linear evolution generated by 
the operator M 2 , which is defined by 

oo 

where 

M nXn = ((e 2 - U + (in + i$(p V )) 2 ) 2 - 3C(*(pi7))) - r, 2 \pd p . 
We want to bound the solution U n of the equation 

e- T d T U ntT = M exp( _ r/2))2)n t/ n)T , (B.15) 
with U n = 1. Recall the definition of L = —p 2 — \pd p , and rewrite -^ e x P (-r/2) 2 n as 

^-p(-t/2) A „ = " (1 + + i*(pe- T/2 )) 2 ) 2 - Xmpe-^ 2 ))^ - e^\pd p 

= (e 2 - (1 + (in + i§(pe- T/2 )) 2 f - %(^(pe~ T /2 ))^j + e~ T p 2 + e~ T L 
= X n (pe- T / 2 )+e- T L, 

where X n (£) = e 2 - (l + (in + i^(0) 2 ) 2 - K($(0) + £ 2 - We want to solve Eq.(B.15): 

e~ T d T U n!T = (e- T L + X n (pe-^ 2 ))U^ T , 

with initial condition U n0 = 1. Observe now that X n is an operator of multiplication by a 
function of prj. Since the commutator [p m , —p 2 — \pd p ] is equal to ^p m ', we find [h(p), L] = 
\ph'(p), and, furthermore, 



e h ^L = (L + \ph'(p))e h ^ 



2 

It follows that the solution of Eq. (B.15) is 



TJ - Je T -l)X n (pe^ /2) T 

as one can check by explicit computation. 

From the explicit form of X n , (in pe 
in Theorem B.l, we see that for any x n e if r , we have 

\\U n , T x n \L r < Cexp(- Co (e^-l)n 4 )e^/ 2 



From the explicit form of X , (in particular, the factor of — n 4 ), and the estimates derived 



\ X n\\q,r 
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Combining this with the Remark of (A.l 1), we immediately obtain 
Lemma B.6. IfU T satisfies 

e~ T d T U T = M exp( _ r/2))2 t/ r , 
with U = 1, then there exist a C(r, q, v) > 0, and a c > such that for any w G H , 

\\U T w\\ q ,r,u < Cexp(-e^/ 2 )|M|^ . (B.16) 

To complete the proof of Theorem 3.4, we also need an estimate of the semigroup 
generated by M v This is simply obtained, however, because from (2.13) we see that 
M v l = M exp (_ T / 2 j 2 j, restricted to functions whose Fourier transform is supported away 
from the origin. Using this fact, and the explicit formula given above for M sxp ^_ T ^ 2 n , we 
see immediately that for any w l £ H q r , there exists a constant c l > 0, such that if U T l is the 
semigroup generated by M v l one has 

ll^^ill,,, < Ce—HH^. 



C. The pseudo center manifold theorem for the singular system 
Eq.(3.8) 

In this section, we prove Theorem 3.7. Before we start with the proof, we wish to point 
out in which sense we are here confronted with a new problem, which does not allow for a 
straightforward application of results from the literature. If we write the system Eq.(3.8) in the 
form 

d T x l = A l x l + N l (x l ,ri,x 2 ,x 3 ) , 
djn = — \ri , 

2 (Cl) 

^t*^2 = ^2*^2 ^2(^1' ^7' *^2' ^3) ' 

d T x 3 = r}~ 2 A 3 v x 3 + T]~ 2 N 3 (x l ,r],x 2 ,x 3 ) , 

then, in view of the spectral properties of Eq.(3.9), there is a "gap" between the "central" part 
(corresponding to x 1 and 77) and the "stable" part (corresponding to x 2 , x 3 ). The problem is that 
we are really dealing with a singular perturbation because the non-linearity in the equation for 
x 3 also diverges as r\ j 0. This problem would be more easily overcome if A 2 were bounded. In 
that case, for sufficiently small 77, the spectra of A 2 and r]~ 2 A 3 would not overlap, and we could 
define first an invariant manifold by "eliminating" x 3 , and then the true invariant manifold by 
eliminating x 2 from the equations obtained after elimination of x 3 . However, since the spectra 
overlap for all values of rj, we resort to a strategy which consists of a converging sequence of 
alternate eliminations of x 2 and x 3 . 
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To define these successive eliminations, we consider two equivalent representations of 
Eq.(3.8), one being Eq.(C.l) above and the other being 

d t x \ = V 2 ( A i x i + N i( x v V,x 2 ,x 3 )) , 

9tV = ~ W ' (C.2) 

d t x 2 = tf{. A 2 x 2 + N l( x V V, x 2i x 3)) > 
dt x 3 = A 3,r] X 3 + ^3\ x l^ V> x 2i X 3J • 

We shall again omit the index rj from A 3 . We obtain Eq.(C.2) from Eq.(C.l) by rescaling the 
evolution parameter of the autonomous system as t + t = exp(r). (Note that time is really 
given by I/77 — t , while we view t and r as the evolution parameters of the vector fields.) 
We will call $^ enter the flow corresponding to Eq.(C.l) and $f ble the flow corresponding to 
Eq.(C.2). A simple inspection of the definition of these flows yields the useful identity: 

K=^(y + t o) ^ x ) = $Z b y^ X ), (C3) 

where 

C = ( x vV), x = (x 2t x 3 ) . (C.4) 

We shall use the relations (C.4) throughout. The identity (C.3) holds for all and for 

rj > 0. Note that the initial conditions are given for the parameter t = and the parameter 
r = log(£ ), and that r](0) = t . Thus, r](0) is small if the parameter t has been chosen 
sufficiently large. (The bounds on the nonlinearities are uniform in £ > as follows from the 
calculations.) 

Let h be a function of ^. This function will always be an approximate invariant manifold 
for one of two problems . To define these problems, we first introduce two effective non-linearities 

F j (h ;£,x 2 ) = N j (x l ,ri,x 2 ,h (€)) , for ,7 = 1,2, 
Gj(h ;^x 3 ) = N j (x v r), h (£),x 3 ) , for j = 1,3. 

We then define two equations (corresponding to the two different time scales Eq.(C.l) and 
Eq.(C2) of the same problem Eq.(3.8)): The first equation will be called the center system: 

d T x l = A l x 1 +F 1 (/i ;^,x 2 ) , 

d T i] = -\r] , (C.5) 
d T x 2 = A 2 x 2 + F 2 {h ; £, x 2 ) . 

Similarly, we define the stable system 

d t x l = r] 2 A l x l +r] 2 G l (h ;£,x 3 ) , 

d tV = -\r? , (C.6) 
®t x T> = A 3 x 3 + G 3 (h ;£,,x 3 ) . 
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Assume now that h 2 and h 3 are two given functions of x 1 and r\. We define a map 

J : 



h 2 \ I h! 2 



h 3 J \h 3 

through the following construction: We let h' 2 (£) be the function whose graph is the invariant 
manifold for the center system Eq.(C.5) with non-linearity F-(h 3 ; £, x 2 ), and similarly we let 
h' 3 (£) be the function whose graph is the invariant manifold for the stable system Eq.(C.6) with 
non-linearity G- (h 2 ; £, x 3 ). Our main result here is 

Proposition C.l. The map J has a fixed point (h 2 , h 3 ). This fixed point provides an invariant 
manifold for the system Eq.(3.8). 

Remark. We shall in fact show that 5F is a contraction in a suitable function space. In 
particular, we show that 9^(0, 0), the n-fold iterate of 5F, converges to the limit (h 2 ,h 3 ). The 
intuitive approach behind this construction is that the 5F n (0, 0) provide a sequence of successive 
approximations to invariant manifolds for the Eqs.(C6) and (C.5), in which the non-linearities at 
the step are given by the approximate solutions for the invariant manifold problem of the other 
equation: The non-linearities are then F^{h 3 ~ l , £, x 2 ) (in Eq.(C.5)) and G^h^ 1-1 ^ x 3 ) (in 
Eq.(C6)). 

Proof. That the systems of equations (C.5) and (C.6) have invariant manifolds follows from 
our estimates (given in Appendix B) on the semi-group generated by the linear operators A 2 
and A 3 , and our estimates on the non-linear terms. (For expositions of this theory that are 
particularly relevant in the present context, see e.g., [H, M, G].) The functions h 2 and h 3 whose 
graphs define the invariant manifolds satisfy well known integral equations, see below. 

Fix h = (h 2 , h 3 ) and consider the Eq.(C.5). We want to find the function h' 2 (h; £) which 
eliminates x 2 . To construct h' 2 , we first consider the equation 

d T x x = A x x x + F l (h 3 ;£,h 2 (£)) , 

a 1 ^ ' 

d T rj = 

This is a differential equation on a finite dimensional space and we let h) denote the 

corresponding flow. (Of course, the ^-component of this problem can be explicitly integrated.) 
We can then formulate the problem of finding the invariant manifold which eliminates x 2 from 
Eq.(C.6) by looking at the map defined by h t— > 3 2 (h) where 

3 2 {h)= f dre- A ^F 2 (h 3 -^-h),h 2 (^-h))). (C.8) 

J —oo 

(A particularly clear derivation of these equations can be found in [G].) In a similar way, we 
define the flow h) for the equation 

d t x x = r ] 2 A 1 x l + rj 2 G 1 (h 2 ; x v h 3 {£)) , 
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and the map 

9 3 (h)= I dte- A ^G 3 {h 2 ^l(^h),h 3 (^-,h))) . (CIO) 

J — oo 

We now specify the function spaces in which we work. Recall that x { G R^, rj G R and that 
£ G R^" 1 " 1 . We let £ c = R^ © R with the usual Euclidean norm. We also assume that £ 2 and 
£ are the Banach spaces in which the x 2 and x 3 live. In our problem, these Banach spaces are 
the Hilbert spaces H r and H qru , but since we believe the present theory of singular vector 
fields may have further applications, we consider the more general case for the moment (see, 
for example, [W2]). These Banach spaces should have the C k extension property [BF]. The 
functions h 2 and h 3 will be Lipshitz functions from a ball of radius r in £ 2 and £ 3 , respectively. 
They satisfy h-{Q) = and are tangent at the origin to £ J , for j = 2, 3. Thus, we define the 
metric spaces, for j = 2, 3: 

Mj, a = {hj-.E'^E* I h j (0)=0,\\h j (t)-h j (i)\\v < *M-i\\}- 
We also define a distance 

\\hA£)-hA£)\\ F 
^ ,..(/',•/',) = V J Ul 3 , (Cll) 

and introduce the notation 

Standard results about the existence and uniqueness of solutions of systems of differential 
equations now imply that 

\\& T fch)-*l&h)\\ < Ce*M\\t-£\\ , (C.12) 

while 

\\^;h)-^;h)\\ < Ce^p^(h,~h), (C.13) 

for any j3 2 > {N — l)/2. Analogous estimates hold for the flow \f/ 3 , though in that case one can 
choose any exponential growth rate (3 3 > 0, provided 1 77 j is sufficiently small. This is due to the 
presence of the factor of r\ A l in the first equation of Eq.(C.6). 

With this in mind we define two more metric spaces (for j = 2,3): 

K m , Dj = {* T : R + x £ c x IK 2)(T x ^ - £ c | 

%(C,h) = Z,y T (0;h) = 0,V T isC 1 inr, 
\\9 T {t,h)-9 T {ih)\\ 

< ^e^ M (iie-lii+P^(^^)iieii)}, 
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with a corresponding Lipshitz metric 

^(*,*) = II*-*!!*,, (c.i4) 

where 

II* I lac,- = SU P SU P 



**(OII 



i>0 ?ee c 



lllll 



These spaces are modeled on those used in [EW]. 



Remark. Since we are interested in local invariant manifolds, we will assume that the non- 
linear terms have been cut off outside a ball of radius r in each of their arguments. Since in the 
applications of this paper all our functions are elements of Hilbert spaces, we can assume that 
there exist smooth cut-off functions which are equal to 1 inside a ball of radius r/2 and are equal 
to zero outside a ball of radius r, and we multiply each of the non-linear terms in Eq.(3.8) by 
such a cutoff. For example, in Eq.(C6), we certainly need to cutoff the function rj 1 by rj 2 x(v) 
(where x is the cutoff function) to avoid blowup problems. 

Given this setup, we show that the map 5F is a contraction of 3"C 2 a x ^-3 a- m terms of 
the notation given above J is now defined as jF(/i) = (3 2 (h), H 3 (h)). One must first show 
that J maps this space to itself. This step is however an easy variant of the argument which 
shows that 5F is a contraction, and we leave it as an exercise to the reader. To show that J is a 
contraction, we use the maps (C.8) and (C.10). Then we see that the "j" component, j = 2, 3, 
of 3(h 2 , h 3 ) (0 - 3(h 2 , ~h 3 ) (£) is given by 

A,- = f_ dre-^fahitrf-UjfctiT)), (C.15) 

where 

U 2 (h;^r) = F 2 (h 3 ;9 2 T (^h),h 2 (9 2 T (^h))) 

= N 2 {^;h),h 2 ^l(C,h)),h 3 {^ 2 T (C,h))) , 
U 3 {h;t,r) = G^^h^i^h))) 

cf. Eqs.(C.5), (C.6). Consider now A 2 . From the estimates on the non-linear term N 2 inEq.(3.8), 
we see that "J 2 is a multi-linear function of its arguments. Thus, we can estimate the difference 
in the integrand of A 2 by the sum of the differences in the arguments of "J 2 , multiplied by the 
Lipshitz constant of J 2 . Because we have cutoff "J 2 outside a ball of radius r, this Lipshitz 
constant can be made arbitrarily small by making r sufficiently small. Thus, calling this Lipshitz 
constant £ 2 (r), we see from the estimates on e Alt which follow from the results of Appendix B 



Nonlinear Stability 



39 



and from Eqs.(C.ll)-(C14) that 

||A 2 || £2 < ^ dr ^e- N ^H 2 {r)[p^{K~h) 

+ \\^ T (t h) - * 2 T (£, h)\\ + \\h 2 (V 2 T (Z, h)) - h 2 (V 2 T (Z, h)) || e 
< j^dr^e- N ^\(r)(p^(h,~h) 

+ p K>r (h,h)Ce* T + p Ma (h,h)Ce^ < const. £ 2 (r)p Ma (h,h) . 

Thus, we have shown that 5F is a contraction. 

We next consider the manifold M given by (£, h 2 (£), h^iC)) — where x l is in a small 
neighborhood of and 77 is in a small positive interval < rj < t] . We want to show that M is 
indeed an invariant manifold for the full system Eq.(3.8). From this it follows, since the flows 
^stable an( j ^center ^ e q U i va i en t^ u p t Q rescaling of time, that M is also an invariant manifold 
for the Eqs.(C.l) and (C.2). If we set x 2 = h 2 (£) and x 3 = h^t;), then the third equation of 
Eq.(C.2) is satisfied because the third equation, when restricted to the manifold x 2 = h 2 is just 
the second equation of the stable system Eq.(C6). with non-linearity G 3 (h 2 ; . . .). To see that 
the remaining equations are satisfied just note that the first, second and fourth equations in the 
full system Eq.(3.8) become, after rescaling of time, 

x x = A^x^ + N l {x 1 , rj, x 2 , :r 3 ) , 

1 3 
V = ~W ' 

and if we set x 2 = h 2 and x 3 = h 3 , we see that we are just on the invariant manifold for 
the center system Eq.(C5). Hence, we have found the invariant manifold for the full system 
Eq.(3.8). 



D. The vanishing of the non-linearity at zero momentum 

In this appendix, we prove Proposition 4.3. This proof is essentially a scaling argument. We 
shall study the nonlinearity N l (x l7 r] 7 x 2 , x 3 ) and we restrict it to the invariant manifold, i.e., 
we replace it by N 1 (x 1 , rf) and let rj go to 0. In particular, we shall show that only one term 
survives, namely the one which is cubic in x 1 , and all others go to as rj — > 0. 

To prove this, we will analyze the nonlinearities N- term by term, using their definitions 
as given in Eqs.(3.6) and (3.8). Recall again that A 1 = since we are considering here the 
projection onto the first eigenvalue of L. In Eq.(3.6), the nonlinearities are given by the terms f 2 , 
f 3 , f 4 , and g, and these have been bounded in Proposition 3.1 and Proposition 3.2. Recall finally 
that every factor of w c contributes a factor of e~ T ^ 2 = rj and every factor of w s contributes a 
factor of e~ T = rj 2 to these bounds. 
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Begin by considering the contribution from f 2 . According to Eq.(A.9), we can extract 
another factor of r\ from Eq.(A.lO), by using the quadratic nature of K 2 , cf. Lemma A.l. 

Using Proposition 3.1 and Proposition 3.2, we see that the only contributions from / 3 , / 4 , 
and g which do not vanish as 77 — > are those of the type (w c ) 3 in f 3 , of the type w s (w c ) 2 in / 4 , 
and of the type (w c ) 2 in g. 

We start by analyzing / 3 . If we write it out, we find 



^ 2 (/ 3 (u> c ))(p) = ^ 2 x( $ (w)) y 

,tT'*(1/2) 

x r? 2 / dp l dp 2 (ptf)®' (p 2 r]) 

i? 7 - 1 $(-l/2) 

x ^(pj w c (p 2 ) «; c (7/ _1 $ _1 ($(p7/) - $(^77) - ^(P 2 »/))) , 
c/ Eq.(2.1 1). Upon taking 77 — > 0, this converges to 

X(0) y dx^ (x)^(x) y dpi^^CPiVW^CP-Pi -P2) • (D- 1 ) 

Analogous arguments can be used to discuss the "surviving" terms of f 4 and g. We just 
summarize the steps analogous to the calculation of f 3 . One gets, as 77 — > 0, 

V~ 2 (f 4 {w c ,w s V,v))(p) -> 6 X (0) y rfa;^ (a;)w e (x)v?o(^) ^ dp'w c {p')w\p -p'; x) , (D.2) 



and 



7/ '(9(w c , w s i],i]))(p) -> -3u e (x)$|(a;) y dp w c (p-p')w c (p') . (D.3) 



We next study these limiting expressions in the basis {^ n (p)}^L °f eigenfunctions of 
L = — jj 2 — 5£»<9 p . Then we can write 7u c (/j) as 



^ (P) + E x ^n(p)- (°- 4 ) 
n=l 



The crucial remark is now that on the invariant manifold, x 2 will be replaced by h 2 n \ and 
similarly w* will be equal to h 3 . We now compute the limiting forms of h 2 and h 3 , and then we 
substitute these values in Eqs.(D.l)-(D.3). Consider the equation for h 3 . Then from Eq.(C9), 
we have 

d t x l = v 2 GM;x lt t4(^)) , 
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because we are considering the case N = 1 where the linear part vanishes. We also have from 
Eq.(C.lO), 

K(x vV )= f cfte-^'Ga^;*?^,*/;^),^^^^;^))). (D.5) 

J — oo 

Now, when 77 = 0, we have (£; /i) = (xj, 0; h) = x v and this reduces to 

^(*i,0)= f dte- A ^G 3 (h;;x v 0,h;(x v 0)) 

J-00 (D.6) 



-00 
1-1 



Note next that for 77 = we have A 3 = M Q , cf. Eq.(3.8), and this means A 3 = Q per L per . We 
denote by £ n (x) the eigenfunctions and by a n the eigenvalues of Q peT L per Using Eq.(1.9) and 
Theorem 1.1, we see that a n = X £=0 n _ l and therefore they are given by a x = —0(e 2 ) and 

a n ?s — (1 — (n — l) 2 ) 2 , when n 7^ 1. Then the n th component (in this basis) of h 3 (at 77 = 0) 
is given by 

(p) = . (_ 3 | dx^KW*)) J dp'w c (p-p')w c (p') , (D.7) 

since all other terms vanish in the limit 77 — ^ 0. We next substitute the value Eq.(D.4) for w c 
and set x 2 = h 2 in Eq.(D.7), and get 

^.(»)(p) = -afr-i . (_ 3 | efe£>K(*)^(*)) 

X (7 dp'V> (p')^o(P-p') + <?(^2 + (^2) 2 )) • 

Next, we replace vf in Eq.(D.2) with h 3 , and in that same equation make the substitution for 
w c that we used above, and we find: 



\%^YjT- x { \ dx£ n (x)u e (x)<pl(x))( dx' (p (x')u e (x')^ n (x')) 

n=0 J J 

x ^ dp, dp 2 i/j (pi)i/j (p 2 )^o(p - Pi - P2) + C(^1^2 + C*: 



(D.8) 



Thus we see that the only terms which survive in N t and N 2 in the limit 77 — > result from 
adding together Eqs.(D.8) and (D.l). We obtain 

X = x\{ j dx<p (x)<pl(x) + 18 / cfa / e B (x / )« e (a; , )rf(x / )) 

^ n=0 ^ 

x(|dx"^> £ (x")^(x"))} (D.9) 
x ^ dp l dp 2 ip (p l ) ip (p 2 ) %(P ~ Pi -P 2 ) 
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This coefficient will turn out to be exactly the same as that which appears below as the coefficient 
of the cubic terms in the center manifold in the periodic case, and since we know that in periodic 
case this coefficient (and indeed, the entire nonlinear term) is zero, it must vanish in the present 
case as well. The only remaining point in the proof of Proposition 4.3 is the computation of the 
coefficient of the cubic term in the equation in the center manifold in the periodic case, and we 
do that in the following subsection. 

Remark. The above argument might seem incomplete since it ignores the 0(x l h 2 + (h 2 ) 2 ) er- 
ror terms in (D. 8). In fact, those terms vanish for x j small. To see why, note that our computations 
of the i] — > limit of f 2 , / 3 , f 4 and g apply also to the nonlinear term N 2 (x v r), h 2 (£), ^(C)) 
in the equation for h 2 in (4.1). Thus, in the 77 — > limit h 2 satisfies: 



Using the estimates on h 2 and derived above, we see that this equation implies h 2 (x 1 , 0) =0 
for all x l sufficiently small, and hence the error terms in (D.8) vanish. 

D.l. The non-linearity in the periodic case 

In this subsection we compute the explicit form of the non-linearity (which we know to be 
because the invariant manifold is made up of fixed points in this case). But this explicit 
form will allow us to compare it with the expression obtained in Eq.(D.9) so that the proof of 
Proposition 4.3 will be complete. 
We start from the equation 



Let y be the component of v in the direction of the highest eigenvalue, cr = 0, of L per , and y, 
the projection onto the directions £ n , defined after Eq.(D.6), associated to the eigenvalues a, 
Then the invariant manifold can be written in the form 



d^*^, 0)^,0) = A 2 h*(x v 0) + N 2 ( Xl Ah* 2 (x v 0),h;(x v 0)) . 



— V 



.3 



(D.10) 



Vn = Y n(Vo) . 



n = 1,2,... . 



(D.ll) 



Using the fact that the eigenfunction with eigenvalue is u' £ , we can decompose v as: 



oc 




(D.12) 



71=1 



the projection of Eq.(D. 10) onto the invariant manifold leads to 



d T V0 



I 



dxu' £ (x) (3u s (x)v(x) 2 + v(x) 3 ) . 



(D.13) 



Note that there is no linear term because a n = 0. 
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We are interested in the exact form of the cubic term in y on the r.h.s. of Eq.(D. 13). There 
are two contributions, one from v 3 , leading to 



-yl J dxu' e {x)\ (D.14) 
and one from the quadratic non-linearity: 

oo „ 

Y = -6y J2Y^(y ) J dxu' e (x)u e (x)u' e (x)i n (x) . (D.15) 



n=l 



Here, Y^ 2 \y Q ) is the quadratic term in y of Y n . Substituting Eq.(D.13) into the equation for 
Y , we find the perturbative result: 



y£ 2) {Vo) = yt-la- 1 J ' dxl n {x)u e {x)u' e {xf . 
Inserting into Eq.(D.15), it is seen to become 

OO „ „ 

Y = -y 3 lSj2^n l J dxu' £ (xfu e (x)Ux) J dx'Ux')u e (x')u' e (x') . (D.16) 



n=l 



Combining Eqs.(D.14) and (D.16), we get the desired result, namely that the cubic non-linearity 
in the periodic case coincides with the quantity X of Eq.(D.9), provided we recall that <p = u' e . 
This completes the proof of Proposition 4.3. 
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